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Abstract 

We outline the relationship between the thermodynamic densities and quasi-particle 
descriptions of spectra of RSOS models with an underlying Bethe equation. We use this to 
prove completeness of states in some cases and then give an algorithm for the construction 
of branching functions of their emergent conformal field theories. Starting from the Bethe 
equations of D n type, we discuss some aspects of the Z n lattice models. 
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1. Introduction 



A number of questions of interest regarding interacting many-body quantum systems 
are posed within the framework of the limit of the number of particles going to infinity. 
Two situations in which will be relevant to the discussion in this paper are the thermo- 
dynamics at non-zero temperatures and the low-energy eigenvalue spectrum of the hamil- 
tonian. To set up the thermodynamic formalism, it is imperative that one has a control 
over counting the states of the system, in order to compute the entropy. On the other 
hand, the low-energy spectrum spectrum is, generically, of a form that is interpreted as 
that describing "quasi-particles," which are then said to populate the allowed energy levels 
subject to composition rules. In order to determine these rules, it is again necessary to 
classify and count the states. 

The main point of this paper is twofold - first, we point out a direct (and innocent- 
looking) connection between the thermodynamic formalism and the quasi-particle descrip- 
tion in a wide class of one-dimensional quantum spin systems whose spectral information 
is expressed in terms of solutions to the Bethe equations. This will then enable us to carry 
out a detailed counting of the states on a finite lattice, and consequently to contruct the 
branching functions of the emergent conformal field theories. 

A quick word about where to place these results: The correspondence between 2- 
dimensional lattice models of statistical mechanics, their 1-dimensional spin-chain coun- 
terparts and their field theory limits - both conformal (massless) and their integrable 
off-critical (massive) counterparts, have long been explored in the context of the particular 
models we consider. In ||, the study of the critical generalized RSOS models was un- 
dertaken using "Bethe ansatz techniques," that had earlier been developed for excitation 
spectra [|Jand for thermodynamics @, (see also 0]@). In addition, the authors of [|J com- 
puted the central charges by taking the T — > limit of the thermodynamic calculation 
of the entropy, following H. These methods have been further extended for the critical 
models of higher rank algebras (J3|, ||13|| ) and for the off-critical models whose masses have 



been determined and S-matrices computed in [El . The key ingredient in all of these stud- 



ies is the densities of the roots of the Bethe equations, and to quote the authors of |T! 
"the set of densities ... determines all macroscopical observables in the model (pg. 307)." 



In [fL6| , the composition rules (flT(J) and the quasi-particle spectrum ( ||14|| ) were 
used to construct the partition function of the low-temperature quantum 3-state Potts 
spin chain, which were shown to to give g-series expansions of known modular branching 
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functions of the associated conformal field theories. The work of |2THused the method of 



20(1 to extract the central charge by taking the q — *■ 1 limit of the g-series, uncovering, in 



the process, the form of the dilogarithm identities very similar to those that had featured 
in the central charge computations via the thermodynamics of these Bethe ansatz solvable 
models. This similarity was exploited to formulate new (/-series identities for branching 
functions in JI7|] , |[8|| , [23 and ]23|] . The work of [ 24 ] served as a bridge between the thermo- 
dynamic calculations of ||, Qand fl3|1and the quasi-particle approach of |l5|]and [p!6| , thus 
"explaining" away the origin of these identities into the structure of the Bethe equations 
themselves. The key step that allowed this connection is the innocuous equation (4.4). 
The proofs of some of these forms were already known (J28|), but since then, proofs of 
some more of these identities have appeared ( [^9[ , [[30| ) . 

A brief outline of this paper: Section 2 explicitly shows the counting procedure that 
was followed in [^0|to get the composition rules for the simplest series of RSOS models. 
Section 3 outlines the procedure for those models based on the simply laced algebras. 
Section 4 sets up the correspondence between the thermodynamic and the quasi-particle 
descriptions. Section 5 outlines the algorithm for setting up the g-series that count the low- 
lying states which are (conjectured, in general, to be) branching functions in the vacuum 
sector of the coset conformal field theory. Section 6 talks about the Z n models of ||31||as an 



example, and the Appendix works out a useful combinatoric identity to count the states 
in the ground state sector of arbitrary spin su{2) models. 



2. Counting states in the ABF models 

The RSOS models that were introduced in Q. In ||, the eigenvalues of the transfer 
matrix of these models were expressed in terms of the roots of the following equation (in 
their notation, I = r — 2): 



sinh (2(]T2)( A J -»)) 



N 



n ^ r 2 sinh(^ T (A J -A fc + 2i)) 



(2.1) 



Qj is a phase factor, iV the size of the chain and / is an integer. 



The solutions to these equations are assumed to be of the "string" form |l|]and follow- 



ing 



we write them in the form 



\ j = \ + i(j + l-2j 1 ), l<ji<j, l<j<l, 3mA = 



0. 



(2.2) 
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These assumptions were made in the bulk in , (with a number of caveats referring to the 
appearance of solutions not of the above form) but here we shall impose these on finite-size 
lattices. (Of course, the same thermodynamic limit can be recovered.) On multiplying out 
the equations for the components of each string, 0§we get a set of equations for the real 
parts, A, of the strings. We then take the logarithm of ( |2.1| )so that the integer branches 
are distinct: 

l M k 

fc=l f=l 

where A^ labels the center of the fi-th string of length j which is a root of eq. ( |2.1|) . The 
functions tj and Qjk are defined below. 

t J (A) = /(A;b'-l| + l), (2.4) 

(min(j,k) — l \ 
/(A; \j - k\) + 2 J2 f(X ] \j-k\ + 2i)+f(X;j + k)j (2.5) 

1 . / sinh lJIT2)( in - X ) 



f(X;n) = In t^Zfi2 1 (2.6) 

M ' ' 2ni Vsinh|^2j(m + A)/ 1 ' 

for integer values of n/(l + 2), and is otherwise. 

For pedagogical reasons, we shall carry out in some detail, the counting procudure 
which is central to all further physical elaborations. We define 

l M k 

z 3 (\) = tjix) - - E e ^( A ' A ")' ( 2 - 7 ) 

k=l v=\ 

so that the (half-) integers ij M satisfy 

Z j (Xi)=27ri I jf. (2.8) 

If we assume that Zj(X) is monotonic, then the range of the integers, Alj = Ij^ max — Ij ; min 
is set by taking the difference of the limiting values, Zj(±oc). A result used repeatedly is 
the following: 

1 / _i_ 9 — r? 

d n = — [/(+oo, n) - f(-oo, n)} = \ , n < I + 1. (2.9) 

Therefore, 

l min(j,k) — l 

A/j = Ndj — ^ Mfc{(l — 5jk)d\j-k\ + dj + k + 2 J2 d \j-k\+2i}, (2.10) 

k=l i=l 
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where min(i,j) picks out the smaller of the two values % and j. Substitution (|2.9|)into 
( p.lOp and performing the sum over the variable i takes the form 



(/ + 2)AI J =N(l + 2-j)-Y, M k {{(l + 2)-\j-k\)- 5 jk {(l + 2) - \j - k\) 



fc=i 



+ ((/ + 2) - (j + k)) + 2((l + 2) - \j - k\) {min(j, k) - l) 
— 2min(j, k)(min(j, k) — l)} 

Splitting up the sum over k to get past the min(j, k) function, we get 



(2.11) 



(/ + 2)Aij = N((l + 2)- j) - J2M k {((l + 2)-j + k) + ((I + 2) - j - k) 

fc=i 

+ 2((/ + 2) - j + k)(k- 1) -2k(k-l)} 

- M 3 {({1 + 2) - 2j) - 2j(j - 1) + 2(1 + 2)(j - 1)} 

i 

- J2 M k {{(l + 2)-k + j)+((l + 2)-j-k) 
h=j+i 

+ 2((l + 2)-k + j)(j-l)-2j(j-l)} 

3-1 

= N((l + 2)- j) - ((I + 2) - j) kM k + 2fM J + (I + 2)M j 



(2.12) 



k=i 



-2(l + 2)jMj-2j ({l + 2)-k)M k . 

k=j+i 

Consider the special case of the strings of length /: setting j = I above, 

i-i 

(/ + 2) Alt =N(l + 2-l)-2(l-l + 2)^2 kM k + (2/ 2 + (/ + 2) - 2(1 + 2)1) Mi 

k=i 

i-i 

= 2N - 4V kM k + Mi(2 - 31), 

tx (2-13) 

= 4(y - J2 kM k) + + M t (2- 30, 
k=i 

= (l + 2)M h 

where we have used the sum rule on the total number of roots in the last step. We thus 
notice the remarkable feature that is enforced by the string hypothesis and the monotonic- 
ity assumptions on Z(X) that All = M;, i.e., all the allowed values for the integers that 



specify the locations of the real parts of strings of length / are always occupied - there are 
"no holes" in this sector || . Going back to (|2.12|) , we replace the Mi occurring in the last 
sum by |r 52fe=j, kM^, and after a few simplifications and dividing by (1 + 2), we get 

1-1 



AI j = N(l-±) + j 1 %2kM k (J-l) + -f J2 M k (k-l) + M^ 

k=l k—j+l 



l-l 



= N(l- 3 1 ) + 2 ^{mm(j, k) - J -^}M k (2.14) 

fc=l 

= N{min(j, 1) - J -^} +2^{mm(j, k) - y }M fc . 

fc=i 

Note that the coefficients in the sum in the last line in ( |2.14j) are the j, kth matrix elements 
of the inverse Cartan matrix of the Ai-i root system! We write the coefficient of N in a 
suggestive way for comparison with later results. 

These integer specifications can be used to label the eigenvalues of the hamiltonian. 
Does this give a complete classification? We shall evaluate the number of ways these 
integers can be ascribed to the A s that determine the eigenvalues, i.e., perform the following 
cominatorial sum: 

{M 3 } j 

To do this we shall follow |7[]. First of all, we reinstate Mi into ( |2.15| )by replacing Mi by 
the sum rule, re-expressing it in the variables M2, . . . , Mi 

/iV/2 + EU(*-2)MA ti /M 3 -+2EU +1 (fc-j)MA 

fa \ w-V^w. )&K m, )• (2 - 16) 

where we notice that M2 now occurs in only two places, and the sum over M2 is of the 
form 

which can be evaluated as the coefficient of x A in (1 + x) c (1 — x 2 ) _B_1 . 
We are now left with the following: 

{M fe },fc>3 J K ' fc=3 v ' 



where the contour is around the origin. The sum over M3 is now of the form X]j>o (^"O^' 
which is just (1 — We can carry on doing this for successive values of k. In 0, 

this iterated sum has been encoded in the following form: 

s 1 U(i±£)^ 1 n 1 (2 . 19) 

2ix V x / 1 — x 2J -- L l — u ■ 

J fc=3 3 

where 



x 
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(uj-1) =u j _ 1 u j+1 , u 3 = ( 1 + x ^!_ x 2y u ? = x ■ ( 2 - 20 ) 

These are essentially Chebyshev polynomials of the second kind, i.e. Uj(x) = Uj(^yP^^), 
where 

sin((n + 1)6) 

U n cos0 = V \ - )V) . 2.21 

SUl</> 

By looking at the poles and residues of Ui(z)/Ui+i(z), we notice that 

UAz) 1 sin 2 (^) 



1 ^ sm (tt^J 



U l+1 (z) / + 2^ z _cos(^ 
Therefore we can express the sum as 

s = * K f JL / i( li£)-^^_ ( , 23) 

1 + 2 2n J \ x J i-(4cos 2 (^)-l)x V ' 

(notice that the square roots cancel out). By deforming the contour to pick up all the 
other poles except at the origin gives us 

2 ^ 2 7Tj_ )(2 ( 7Tj_ N 

1 + 2^ K l + 2 JK y l + 2 !! ' v ! 

i=i 

which is precisely the multiplicity of the singlet in the iV-fold tensor product of the funda- 
mental representation of U q (su(2)) at q = exp(^ii). This is the expected number because 
of the the equivalence |38|]of the construction of the model state space to truncated tensor 



product representations of quantum groups at roots of unity. (This method of performing 
the sum, which is lifted from 0, where it was used for the I — > 00 case, was deployeded 
for the I = 4 case in |H|. See also [j39|.) 
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3. For the simply-laced algebras. 

A bit of history: In H, higher spin representations of su(2) were studied. From the 



su(2) models, we can generalize even further, by considering the models of |TT[as solved in 
0, by starting with the corresponding Bethe equations. In fact, the authors of 0made the 
observation that the structure of the Bethe equations could be cast in a way that could be 
generalized to the case of all simply laced algebras even though the models whose spectrum 
these equations would parametrize were not known. This game was further extended to 



the elliptic case in ||12|| and by Kuniba ||13|| who set up the thermodynamics of (hypothetical) 



systems whose energy eigenvalues were parametrised by solutions to the (trigonometric) 
Bethe equations associated with all (i.e. not just simply-laced) untwisted affine algebras. 

We shall briefly outline the counting procedure for the simply-laced cases which closely 
follows section 1. This will highlight the key feature characteristic of the connection be- 
tween the quasi-particle description (which is a recasting of the map between the integers 
and momenta in a suggestive language that posits this physical paradigm) and the ther- 
modynamic formalism. 

For the simply-laced algebras, the Bethe equations are of the form 0: 



sinh(^(A^ + is6 ap )) j N _ fl(a) * ^ sinh(£(A< a) ~ ^ + ^C ah )) 
sinh(^(A( a) -z^ ap ))J 3 .-^ v*' 



,(«) 

JJ^sinh^Af -\^-iC ab) ) 

N a = Ns[C-i] ap , (3.2) 

where iV is the size of the lattice C is the Cartan matrix, L = I + g, for integer I, g is the 
dual Coxeter number and s characterizes the type of fusion. In what follows, the string 
hypothesis (as in ||) can be written as: 

\f =x + i(j + l- 2ji), l<ji< j, (3.3) 

A is real and denotes the center of the string and j denotes its length, which is further 
assumed to satisfy 1 < j < I. We shall therefore impose Y^j=i = N a , where 

denotes the number of strings of colour (a) , length j and N a is given by (|3.2|) . As before, 
we multiply out the Bethe equations for the components of each string, and end up with 
equations for the real parts of the roots. We then take the logarithm of the multiplied out 
fl3J])so that the integer branches are distinct: 



Ntft&M) = 2*il% e$\$ a) ~ X t ib) ), (3-4) 

6=1 fc=l u=l 
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where \j} labels the center of the /i-th string of length j and color (a) which is a root of 
eq. (|3.1|). The functions t* a s and are defined below. 

min(j,s) 

t^ s W = S ap J2 m\j-s\+2k-l), (3.5) 

k=l 

min(j,k) — l 

Jab), 



e)T (A) = s ab ( /(A; \j-k\) + 2 J2 \i - k \ + 2z ) + /( A ^ J + k ) 

1=1 

in(j,k) 

-I ab [ /(A;|j-*|+2i-l) 



(3.6) 



i=i 

/<A:»)= — iu( •;t?!/V , (•>•-) 



1 /sinh i^fm — X) 

— In 

27rz \sinh||(m + A) 



for integer values of n/L, and is otherwise. I ab is the incidence matrix of the respective 
Dynkin diagrams. As in section 2, we define the corresponding Z(\), this time with an 
extra colour index. Assuming monotonicity, we find that for strings of length /, the range 
of integers AI^ coincides with the total number of strings of that length, M[ , for all 
colour labels (a). There are no holes in the distribution of integers for these strings. Since 
the Z-strings do not contribute to the counting of states, we can eliminate Mj using the 
sum rule (|3.2|). The ranges of integers associated with the other strings are: 



AI^ = NS^J.^Mf -^(C 6 ).KJ,Mf, j < I, (3-8) 

6=1 k=l 

where Cg is the Cartan matrix of the (simply-laced) Lie algebra Q and (as before) i is 
the inverse Cartan matrix of Ai-±. Note that this specializes to the case of su{2) considered 
earlier, for 8 = 1. This gives a combinatorial count of the number of states of the form 
described above, 

/A/ (a \ 

£ n U, • <™> 

We shall use the equality of binomial coefficients 

b)-(/-s) 

to define a new variable 

N (a) = A/ (a) _ M W ( 3 _ n ) 



which counts the number of holes. In these variables, the allowed range of integers takes 
the form: 

6=1 fc=l 

The observation that AI^ = had been first made in the context of thermo- 

dynamic densities studied in |J, ||. The striking difference is, of course, that in the 
thermodynamic formalism, N is a large number, and only a single order in N is retained, 
and there was no reason to expect this to be true at any finite N. (That the authors 
of [|| notice deviations from the string hypothesis ( |3.ij| ) makes it all the more surprising.) 
We shall exploit this coincidence in order to do away with the complication of choosing 
branches as in ( |3.4| ). 



4. The thermodynamic formalism. 

The procedure of setting up the thermodynamics has now become standard, going 
back to 0. We shall reproduce the equations in ||13|| where the thermodynamics of systems 
based on the Bethe equations corresponding to all simply-laced algebras were studied. 



sinh(^(Aj a) + i(su p \a a ))) 



-i N 



oMfTfT sinh (ir(^ a) -Ai b) + zKK))) 

"j 1111 i / 7r /,(a) ,(a) , N xx • 



sinh (^( A i - i(su p \a a ))) \ b=1 k=1 sinh(^(A} aj - Ajf - i(a a \a b ))) 

N^NsiC- 1 }^ (4.2) 

where N is the size of the lattice a a are the simple roots, uj p are the fundamental weights, 
(■|-) is the canonical bilinear form on the dual to the Cartan subalgebra, C is the Cartan 
matrix, L = I + g, for integer I, g is the dual Coxeter number and s characterizes the type 
of fusion. 

The string hypothesis for the solutions of ( |4.1|) is 

\^ =\ + it~\j + 1 - 2ii), 1 < h < 3 (4-3) 

with t~ x = (a a \a a )/2. A is real and denotes the center of the string and j denotes 
its length, which is further assumed to satisfy 1 < j < t a l. We shall therefore impose 
^2j=i jM^ = N a , where denotes the number of strings of colour (a), length j and 

iV a is given by (O). 



9 



To make the idea of upgrading bulk quantities to finite lattices concrete, we define 
densities for strings and holes ||by requiring their integrals over all A to be exactly equal 
to the number of strings and holes divided by N, the size of the system. This is the precise 
statement that rules the interpolation between the bulk and finite lattices. In the notation 

of m, 

pf (0) = Mj a) /N, of ] (0) = Nf /N, (4.4) 
where /(0) denotes the Fourier transform of / for zero argument: 

/OO /"OO 
f(X)e-^d\, /(0) = / f(X)d\. (4.5) 

One would expect this correspondence to be only true modulo correction terms vanishingly 
small compared to N. In the thermodynamic limit, because solving the integral equations 
involves taking derivatives on the densities, one can choose the branches of the logarithms 
(of the Bethe equations) freely. What we see is that there exists a choice of branch on 
the finite lattice that makes the translation of statements true in the bulk to finite lattices 
possible. 

The equations governing the densities that determine the thermodynamics of these 



systems as studied in [|i~3||are 

r t b l-l 



Mg' m = ^ + EE M ab AT k p { k\ (4-6) 

6=1 fc=l 

where 



£{j)mk _ UJ)km _ smh(min(m/t a , k/t b )x) sinh((J - max(m/t a , k/t b ))x) 

A ab — A ba — 



smh(x/t a i,) sinh( Jx) (4 7) 



Mab = M ba = ^C ab + 25 ab (cosh(x/t a ) - 1), 

tab 

and t ab = max(t a ,t b ). 

Setting x (the Fourier transform variable) to zero, and applying the definition ( |4.4f) to 
(0), we get 

N$ a) = N5 ap [C- A ] pi J JS -EE^f- (4-8) 

6=1 k=l 

where 

K ab = (a a \a b ){mm(t b j,t a k) - —}. (4.9) 
10 



We can also solve for M- in terms of to get 

=NJ2 {K-^[C- A ] pi J k9 - E EV 1 )^, (4.10) 

k=l 6=1 k=l 



where K 1 is defined by 



tj-1 



E E COS^S* = S abS jk . (4.11) 



c=l m=l 



In the above, the index t a l never shows up because there are never any holes in this sector, 
and the sum rule (|4.2| )is used to eliminate the dependence on Mr}. Note that for the 



L t a l 

simply-laced cases, K factorizes into a "level" and a "rank" piece, and ( fOf ) reduces to 
( |3.8|) using ( 3.11| ). The range of integers for the finite system, under this correspondence, 
is A/] a) =Mj a) +iv] a) . 



5. Quasi-particles and (/-series. 

Now for some physics. The first thing to remember is that the hamiltonians are 
connected to transfer matrices of two dimensional statistical mechanical systems with 
positive-definite Boltzmann weights, and for either of the ± signs the ground state has to 
be unique by the Perron-Frobenius theorem. Our classification (assumed complete) via 
the partition of the roots of ( |4.1| )into stringy clusters for the states in this sector (where 
the ground state is expected to lie) indicates that the only states that meet this criterion 
are given by the partition for which the summand in binomial coefficients gives 1. Using 
the variable variable which counts the number of holes in the allowed distribution of 
integers, and look for solutions of = 0. The observation that the range of allowed 
integers for strings of length t a l was equal to the number Mt a i itself indicates that the 
state in which all the roots are of this type is one such candidate for a ground state. The 
only other state is then the state which has a distribution of strings given by 

t p i-i 

Mf =Nj2(K- 1 )H[Cjl_ 1 } k8 . (5.1) 
k=i 

We thus have the two ground states for the two hamiltonians that differ by an overall sign. 
Any other partition of the roots then corresponds to an excited state in either model. 
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A/] a) = N5 ap [C- A ] ] js + - £ $5 ^Mf , (5.2) 



and 



6=1 fc=l 



= at x: (josp^j*- + at]-) - e eV- 1 )^ 6 } ^- 3 ) 

fe=l 6=1 k=l 

encode the taxonomy of all the states of the system except the ground state which has 
been "subtracted out" We can then regard (|5.2| )and (|5.3| )as the fundamental equations 
that determine the excitation spectra of the two models, with the numbers M^' and 
counting the number of these excitations. 

A generically striking property that the spectrum of a many-body quantum system 
(such as the ones described here) possesses is one described as a quasi-particle form. This 
refers to the additive decomposition of the eigenvalues thus: 

lim E-E GS = E e >i)' ( 5 - 4 ) 

a,rules j 

where E refers to the eigenvalues of the hamiltonian (Eqs being the smallest eigenvalue) 
and N refers to the size of the system. The functions in the summand are the dispersion 
curves of what are defined to be quasi-particles of different species a, and the index j runs 
over their number n a . The momentum of a particular eigenstate is 

p - p GS= E I>i- ( 5 - 5 ) 

a,rules j 

The "rules" determine the composition of the eigenstate in question, and in particular, 
could depend on the statistics of the excitations. 

The n a quasi-particles of ( |5.4| )and ( J5.5| )are the excitations which number and 
above the two vacua in the thermodynamic limit. The completeness argument on the 
finite lattice (which we have demonstrated only for a subset of the cases under considera- 
tion) indicates the completeness of the particle picture of the (separable) Hilbert space of 
the emergent field theory. 

The momentum eigenvalue ( |5.5p is written as a sum of the logarithm of the term on 
the left hand sides of ( f4.1| ). Thus, the sum of the integers for any set of roots gives the 
total momentum (with the appropriate taken into account) of the state that these roots 
correspond to. Therefore we see that for every term in AI^ = M- + for which 
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the term proportional to N does not have support, the contribution of these strings (or 
holes) to the momentum of the eigenstate at order one (N ) must necessarily vanish. For a 
system with no mass gap the energy contribution must consequently be zero to order one. 
However, they could (and do) contribute to the spectrum at order 1/N. These have been 
termed "ghost" excitations in [j!6| . The strings (and holes) for which AI^ ~ N have an 
extensive single-excitation Hilbert space dimension and constitute the order one excitation 
spectrum, and the coefficient of N encodes the Brillouin zone scheme of these particles. 



(The relationship of these spectra to their fractional statistics |25|has been mentioned in 
||15|| and in ||26|| , where the realtionship of the statistical exclusion coefficient to the central 
charge of the conformal field theory has been shown.) 

It is seen in numerical studies (which provided the background for ref. ]T^] ) , that 
within each class of states with a particular root content, the sum of the absolute values of 
the integers gives a good estimate of the (approximate) degeneracy of the levels, i.e., those 
states with the same value for the sum of the absolute value of the integers had energy 
eigenvalues that were almost equal. This is reminiscent of the conformal field theory 
definitions of energy and momentum being the sum and differences of the eigenvalues of 
Lq and Lq. One could then presume that this degeneracy would become exact in the 
thermodynamic limit, and that these integers keep track of, in a rather robust fashion and 
on a finite lattice, a classification of states that take on significance (in terms of dynamical 
symmetries) only in the thermodynamic limit. We can then hope to count degenerate 
states in the iV — > oo limit by keeping track of these integers. (For more discussion on 
quasi-particles and state counting in the thermodynamic limit, see ||19|| .) 

In other words, the binomial counting may be refined by introducing a grading, i.e., 
by introducing a variable q, whose powers are the integers. Consider, for example, any 
3 of the 5 integers in the set {±2, ±1,0} can be chosen such that their sum takes values 
in the set {±3, ±2, ±1, 0}. The multiplicity of their occurence can be observed from the 
coefficient of the appropriate power of q in 



jjjj =q~ 3 + q~ 2 + 2q~ 1 + 2 + 2q + q 2 + q 3 



where 

[A + B 



B 



B 



(5.6) 



(5.7) 
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are defined in terms of g-binomial coefficients, 



(q; o) 



A 



(q;q)A-B(q;q)B : 



(5.8) 



where 



(q;qU = i[(i-q J ), 

(non-zero only for integers, A and B, with < A < B). 
Construct the following: 



(5.9) 



fp= e n ^"" Mj<0)7Vj(0)(Mj<0)) 

{M (a) }i.( a ) 



Mj o) + ivj a) (Mj o) )- 



M 



(a) 



{7V] a) }i:(«) 



(a) 



(5.10) 



where the parantheses are used to indicate that we write / p only in the Mj 8 ' variables 
and //j only in the A^- variables. Since large values of A correspond to small momentum 
contributions and therefore, low energy states, and the integers have a monotonic depen- 
dence on A, the largest integers correspond to the lowest energy. We therefore assign zero 
momentum to the term in AI^ proportional to iV (which is the edge of the Brillouin 
zone) in the power of q in the factor pre-multiplying the (/-binomial coefficients. 

We can now count the low energy states in the thermodynamic limit, iV — > 00. For 
the massless case, this is equivalent to computing the partition function of the theory in 
the particular sector, with q = exp(-j^), in the limit iV — > 00, T — > 0, NT finite, v being 



14 



the speed of sound. In the iV — > oo limit, the /„ and fh are re-defined as follows 



Xv 



lim f p 



lim 

N—>oo 



E 



q -±M-K-M x 



X 



n 



J, « 



2 u ° 



[js-TMj ' IV ab lvl k 



M 



(a) 



Xh 



Jim fh 



lim 

N—>oo 



E 



X 



X 



n 



3, a 



N 



(a) 



with the condition that 



lim 

N— >oo 



'iV' 
m 



(?;?)t 



(5.11) 



(5.12) 



These are conjectured to be the branching functions (in the vacuum sector) corre- 
sponding to the conformal field theories constructed as cosets (a la GKO) f2"T|j . Some of 
these have been conjectured earlier (see [Jl^jfor references). The branching functions these 
(/-series expansions correspond to are listed in ||13|| . The central charges which can be ob- 
tained by taking the q — > 1 limit as in |20|| . ||21|| and |18[are also listed in [Oj. Thre exists a 



proof of the su(2) formulas in |2^]for the M-type variables (x P ) and in []30[]for the iV-type 
variables (xh)- 



6. The Z n models and some "Snake-Oil". 



In JT^] it was pointed out that the g-series corresponding to what in our notation would 
be Q = D n , s = 1 and p = 1 and in the M variable coincided with the vacuum character 
of Z n parafermions (with the colour index (a) taking values upto the rank, n). 

We have obtained [^4|this g-series from a Bethe equation, without any reference to a 



physical model. The Fateev-Zamolodchikov model |31]in its scaling limit is known to give 
rise to the parafermionic conformal field theory. The Z3 case was studied in ||10|| [pa using 
the formulation of and we saw that the formulas in section 1 reproduced (modulo 
an orbifolding which we shall indicate later) those results (for the Q = sector of the 
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3-state Potts model). To confirm whether this counting procedure in terms of the roots of 
(pTT|) gives the total number of Q = states, we need to execute a corresponding bimomial 
sum. 

We shall number the nodes of the Dynkin diagram of D n such that the point of 
bifurcation is at 3 and the ends of the fork are at 1 and 2. The expression we are going to 
evaluate is 

_ /M + m 3 \ / m 3 \ fmi + m 2 + m 4 \ /m 3 + m 5 \ /m n _ 2 + m n \ (m n -\ 
M ~ \ 2m! J \2mJ \ 2m 3 J \ 2m 4 J " ' \ 2m n _i ) \ 2m n 

Recall that for the critical 3-state Potts case, which is the n = 3 case of the Z n models 
of Fateev and Zamolodchikov, the number of genuine 0(1) quasi-particle excitations were 
both even and odd in number unlike the corresponding orbifold-related RSOS model which 
only had an even number in keeping with the Z 2 invariance. The even sector only picks 
up the states with the Z 2 positive charge states. We therefore expect that to recover the 
complete Z n charge zero sector we will have to sum over mi half integral as well, which 
will necessarily make the associated m 2 to be half-inegral too. 

Following the elegant and extremely useful method advocated by H. Wilf called 



the "Snake Oil" method ||32| , we define 

/(x) = ^x M S M , (6.1) 



M 

interchange the order of summation and perform the sum over M first. We shall follow the 
convention that the binomial coefficient („) vanishes if B < or if A is a positive integer 
less than B. Using the formula (valid for \x\ < 1), 

we get 

. ^ M /M + m 3 \ x 2m 1+a -m 3 

2-^ X V2mi + oJ ~~ (1 - x) 2mi+a+1 ' 

M>0 v 1 7 v 1 

(Here a = 0, 1 which we have introduced in order to keep track of the odd/evenness of 
the number of excitations.) Plugging this result into the sum over mi, we then need to 
perform 

x -m 3 +a x 2 x mi (m\ + m 2 + m 4 + ct N 



(1 -x) a + 



Ell X \2\ m l ( m l + m 2 + m 4 + OA 
^T^ } > V 2m 3 J 
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which is again of the form ( |6.2| ). After some simplification, the sum over m,2 then takes 
the form 

( x ' 1 - 2x (1 - 2x) 2m 3 2^ x > \2m 2 + a)' { ' } 

The presence of the as alongside the summed (dummy) variable indicates that the 
sum is performed over all (even and odd) the integers, and thus we have 



, , s _ >T-^ ^ 1 - X , X x 2m 3 , 1 - X , 2m 4 fm 3 + m 5 \ 

n x > ~ l^ - l^ i-2x'l-2x^ ' x ' V 2m 4 j 

m 4 + m 6 \ /m n _ 2 + m n \ /m n _i 
2m 5 / " ' V 2m n _i / V 2m n 



x 



(6.4) 



Define fs = (x/(l — 2x)) , so that upon performing the sum over 7713, 
J[X) ~ 1-2x1- f 3 2s h [ x l-h> V 2m 5 

^ J m„,...,m 4 ^ 

m 5 + m 7 \ fm n -2 + m n \ (m n -\ 
2m 6 J " \ 2m n _i / V 2m n 



x 



(6.5) 



Defining 

X 1 — J 3 1 — 6X 

and iterating this process, we end up with 

1 - r -A- 1 

f(x) = — — — TT — —r, (6.7) 
JK ' l-2x 11 1- U K 1 

3=3 JJ 

with 

(J/ 1 - I) 2 = /Afi+l. ( 6 - 8 ) 

and the initial two values of fj for j = 3, 4 are defined above. It can be checked that a 
solution to these recurrence relations with the given initial data is given by 

Plugging it all into ( |6.7|) , we get 

fix) = IzJ^zDE. = J2n M -'x M , (6.10) 
1 — nx *s 

which precisely counts the number of states in the (Z n charge) Q = sector of the Fateev- 
Zamolodchikov spin chain. 
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Going back to the characters, we now consider those in the N variables, which we now 
expect are those that describe the conformal field theory of the anti-ferromagnetic chain. 
In flSfl and |2^1 , it was conjectured (based on expansions of the formulas (|5.11|) as (/-series, 



using M athematica™ ) that these characters are those at radii y/^ on the Gaussian line 
of the moduli space of c = 1 theories. This expectation has been borne out in [|34|]for odd 
n. 

7. Discussion. 

In order to describe the physics that the solutions to the Bethe equations encoded, 
we invoked the positivity of the Boltzmann weights of non-existent models. For the 
anti-ferromagnetic D n models, the classification scheme seemed to describe the states of 
the Z n models of ]3l|] whose Boltzmann weights were, however, not positive for the anti- 
ferromagnetic regime (of the associated spin chain). This condition of "confinement of 
holes" served to rearrange the classification of states in a useful way giving correct results. 
(In passing, let us note that in E3] , the E 8 structure of the Bethe equations emerges owing 



to a similar occurence of a " 'frozen' Dirac sea.") 

The more interesting questions that remain to be answered involve cases where the 
dependence of the energy functional on the roots of the Bethe equations is non-trivial, 



as in the case of the superintegrable chiral Potts model [36], where similar completeness 



studies have been made [37 . 
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8. Appendix 

In this appendix, we shall evaluate the sum over the products of binomial coefficients 
for the case of su{2), arbitrary spin s, by the "Snake-Oil" method. We shall warm up to 
it by first re-doing the sum for the s = 1 case. Let us write the combinatorial summand in 
terms of the "hole" variables (we shall change notation, replacing Nj by rrij, and N/2 = 
M): 

s _ \ - (M + m 2 \ /mi + m 3 \ (m 2 + mA /m n _ 2 + m n \ /m n _A 
M ~ ^ \ 2 mi ){ 2m 2 J\ 2m 3 2m n _i J V 2m n J ' 1 ' J 

mi,m2,...,m„ x /v / \ / \ / \ / 

Next, define 

/(.) = ^x M S M , (8.2) 

M 

interchange the order of summation and perform the sum over M first: 



M (M + m 2 \ _ ^ x _ m2xM+m2 (M + m 2 



v x M 2 = y 

^ \ 2m 1 J ^ V 2mi 

M>0 v 7 M>0 v 

We shall follow the convention that the binomial coefficient (^) vanishes if B < or if A 
is a positive integer less than B. Using the formula (valid for |x| < 1, 

y PV = r JL -r(— ) B . (8-3) 

we get 



x - a 7712 mi /mi + m 3 \ (m 2 + mA /m n _ 2 + m n \ /m n _i\ 
/W ~ (l-^)" 1 V 2m 2 A 2m 3 J" A 2m n _ x ){ 2m J ' 



mi ,?ri2 ,...,m 



2 

where we have defined yi = (cc/(l — x)) .In the same way, we can now perform the sum 
over mi and be left with 



x - a m2 j/!"" 3 , yi ^ 2 /m 2 + m 4 \ /m n _ 2 + m n \ /m„_i 
m 9 m V 2m 3 ;-'A 2m n _! J \ 2m n 

1 2/i~ m3 m 2 /' m 2+TO4\ f m n - 2 + m n \ f m n -! 



^ a-xHi-uir 

m 2 ,...,m„ 

where 



7 / 

(1 — x) (1 — y\) 2 V 2m 3 7 V 2m n _i J \ 2m 



n 



(8.4) 



y^i^rjv^) 2 - ( 8 - 5 ) 

v v^(i -yi) 
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The pattern is thus evident and the process is repeated until we end up with 



1 

where 

V-i = 1,2/0 = x, and (yj 1 - l) 2 = y-\y~+ v (8.7) 
These are the same recurrence relations as before, only this time, 

-1/2 TT ( 1 x 



2v^ 



At this point it is easy to see that if we were to consider any other representation, 
not necessarily s = 1, there would be a term involving M in the binomial coefficient at a 

th 

different place, the s place. That is there would be a factor of the form 

m s -i + m s+ i +Af 
2m s 

(where Af is really equal to N/2), which would give rise to a power (y s _i) _A ^ picked up 
while redifining the summation variable so that the required sum is simply the constant 
term (coefficient of x°) in 

/(*■) = t/ l( J, )u^( J, )- f/ " +l( * ) 



2V3F' ' v 2V5'E/„+2(jl-) 



Once again, using (|2.22|) and closing the contours over the zeroes of the denominator 
we arrive at 

Sn = ^—Ysin\p-)U^ (cos(-^-)), (8.9) 
1 + 2 ^ K l + 2 J 3 v K l + 2 JJ ' K ' 

3=1 

with i = n + 1. This specializes correctly to the s = 1 case evaluated above. 

As is immediately obvious, this method can be used for evaluating the sums for cases 
where the quadratic form Kfj involves only the i,j or the a, b indices. The other cases 
require the introduction of multiple variable generating functions and multidimensional 
residues are required. For more results and techniques associated with counting solutions 



of Bethe equations, see [^Of and pi ]. 
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